Transverse depinning and melting of a moving vortex lattice in driven periodic 

Josephson junction arrays 
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We study the effect of thermal fluctuations in a vortex lattice driven in the periodic pinning of a 
Josephson junction array. The phase diagram current (7) vs. temperature (T) is studied. Above 
the critical current I C (T) we find a moving vortex lattice (MVL) with anisotropic Bragg peaks. For 
large currents I ^> I C (T), there is a melting transition of the MVL at Tm(I)- When applying a 
small transverse current to the MVL, there is negligible dissipation at low T. We find an onset of 
transverse vortex motion at a transverse depinning temperature T tr (I) < Tm(I). 



PACS numbers: 74.50+r, 74.60. Ge, 74.60.Ec 

The interplay between the periodicity of vortex lat- 
tices (VL) and periodic pinning potentials in supercon- 
ductors raises many interesting questions both in equi- 
librium and in driven systems 0. Experimentally, 
periodic pinning has been realized in artificially fabri- 
cated systems like: superconducting wire networks |Q, 
Josephson junction arrays thickness modulated su- 
perconducting films p|, magnetic dot arrays Q and sub- 
micron hole lattices H in superconductors, commensu- 
rability effects in the ground state vortex configurations 
lead to enhanced critical currents and resistance minima 
for "fractional" and for "matching" (i.e. commensurate) 
vortex densities where the VL is strongly pinned. Un- 
der the effect of thermal fluctuations, it is possible to 
have a depinning phase transition of these commensu- 
rate ground states at a temperature T p and a later melt- 
ing transition of the VL at a temperature Tm ||- For 
high vortex densities (i.e. strongly interacting VL) both 
transitions coincide, T p — Tm, while for low vortex den- 
sities both transitions are different with T p < Tm- Out 
of equilibrium, many recent studies have concentrated in 
the related problem of the driven VL in the presence of 
random pinning Js — 1 3 1 . The nature of the fastly moving 
vortex structure for large driving forces has been under 
active discussion lately [jB-JlO motivating both experi- 
mental [jll| and numerical |b|-[l4|] studies. In particular, 
the interesting concept of transverse critical current has 
been introduced in Ref. P]. After applying a current in 
the direction perpendicular to the drive, a finite trans- 
verse critical current I c j r may exist, at least at T = 
P|. For T > 0, there is a very small but finite transverse 
linear response JlO[ ] , with a possible sharp non-linear volt- 
age increase at an "effective" I c ^t r P,|lO|]. For periodic 
pinning, the physics of the driven VL has been studied 
numerically only at T — 0, where a complex variety of 
dynamic phases has been reported @. In this case, it is 
clear that a finite I c ^ r will exist due to commensurability 
effects |l(|, and it has been obtained in jjj for T = 0. 

In this Letter we study the effect of thermal fluctua- 
tions in a driven VL in a periodic pinning potential. In 
our case, the periodic pinning is provided by a two di- 



mensional Josephson junction array (JJA) p|, p^|Jl^ ]. We 
obtain a phase diagram as a function of the driving cur- 
rent (J) and temperature (T), which is shown in Fig. 1. 
For low currents, we find that the depinning and melt- 
ing transitions are separated with T P (I) < Tm(I)- More 
interestingly, for large currents we find an analogous se- 
quence of transitions but for the transverse response of a 
fastly moving VL. Wc find that there is a novel transverse 
depinning temperature T tr below the melting transition 
of the moving VL, T tr (I) < T M (I). 

The current flowing in the junction between two super- 
conducting islands in a JJA is modeled as the sum of the 
Josephson supercurrent and the normal current |15|]l(| : 



J^(n) = 7 sin ^(n) 



$o 00 M (n) 
2itcRn dt 



+ ^(n,t) (1) 



where Iq is the critical current of the junction between 
the sites n and n + /i in a square lattice [n = (n x , n y ), 
fx = x, y] , Rpf is the normal state resistance and #^(n) = 
0(n + /i) - 0(n) - A^n) = A (tl 0(n) - A^fn) is the gauge 
invariant phase difference with A^(n) = ^ jj^ + ^ a A ■ 
dl. The thermal noise fluctuations rj^ have correlations 
(»fc(n,*)iv(n / ,t / )) = S^Sn^Sit - if). In the pres- 
ence of an external magnetic field H we have A M x 
A^n) = A x {n) - A x (n + y) + A y (n + x) - A y {n) = 2tt/, 
/ = Ha 2 /$ and a is the array lattice spacing. We take 
periodic boundary conditions (p.b.c) in both directions in 
the presence of an external current I ext in the y-direction 
in arrays with L x L junctions Jl4| . The vector poten- 
tial is taken as A^n, t) = A a ^(n) — a^t) where in the 
Landau gauge A°(n) = —2nfn y , Ay(n) = and a M (t) 
will allow for total voltage fluctuations. With this gauge 
the p.b.c. for the phases are: 6(n x + L,n y ) = 0(n x ,n y ) 
and 0(n Xl n y + L) = 9(n x ,n y ) — 2nfLn x . The condi- 
tion of a current flowing in the y- direction: ju(n) = 
IextL^S^.y determines the dynamics of a p (t) |TJ]. Af- 
ter considering conservation of current, A M • /^(n) = 
7 M (n) - I^(n - (i) = 0, we obtain: 



, 80(n) 



= -A^-[5 M (n)+r ?/1 (n, t)] 



(2) 



1 



^ = W MlW - ^ $>„(n) + ^(n, t)] (3) 

where 5^(n) = sin[A M 0(n) — A°(n) — we have nor- 
malized currents by Jo, time by tj = 2t:cRnIo/^q, tem- 
perature by Io&o/2TtkB, and the discrete laplacian is 
A20(n) = 0(n+x) + 6>(n-x) + 6>(n+y) + 6>(n-y)-46>(n). 

The Langevin dynamical equations (2-3) are solved 
with a second order Runge-Kutta-Helfand-Greenside al- 
gorithm with time step At = O.lrj and integration time 
IOOOOtj after a transient of 5000rj. The discrete lapla- 
cian is inverted with a fast Fourier + tridiagonalization 
algorithm as in jj^ j. We study the following proper- 
ties: (i) Superconducting coherence: we calculate the 
helicity modulus in the direction transverse to the cur- 
rent = £ (£ n cos 6 X (n)) - ^{(E n sin^(n)] 2 ) - 
(E n s ' n ^( n )]) 2 }' P n or< ier to calculate the helicity 
modulus along x, we enforce strict periodicity in 9 by 
fixing a x (t) — 0] . (ii) Transport: we calculate the time 
average of the total voltage V = (%(*)) = (da y (t)/dt) 
(voltages are normalized by R n In)- (hi) Vortex structure: 
we obtain the vorticity at the plaquette n (associated to 
the site n) as 6(n) = — A M x nint[0 Al (n)/27r] with nint[x] 
the nearest integer of x. We calculate the average vortex 
structure factor as 5(k) = (\-rz J2h K**) exp(ik • n)| 2 ). 

We study JJA with a magnetic field corresponding to 
/ = 1/25 and system sizes of L x L junctions, with 
L = 50, 100. The ground state vortex configuration for 
/ = 1/25 is a tilted square- like vortex lattice (VL) com- 
mensurate with the underlying periodic pinning poten- 
tial of the square JJA (see M])- The structure factor 
5(k) has correspondingly delta-like Bragg peaks. For 
this value of / we find an equilibrium phase transition 
at T c fss 0.050 ± 0.003, which corresponds to a simulta- 
neous VL depinning (corresponding to the onset of resis- 
tivity and vanishing of helicity modulus) and VL melt- 
ing (corresponding to the vanishing of Bragg peaks); i.e. 

T c = T P = Tm- 

First, we have calculated the current-voltage (IV) char- 
acteristics for different temperatures. At T = there is 
a critical current of I c (0) = 0.114 ± 0.002, which corre- 
sponds to the single vortex depinning current in square 
JJA [l8| ]. Above I c (0) there is an almost linear increase 
of voltage until / w 1 where there is a sharp rise of V be- 
cause all the junctions become normal. Similar behavior 
has been reported for T — IV curves for low values of / 
fL9f . We restrict our analysis for currents I < 0.4, where 
the collective behavior of the VL is the dominant physics. 
For temperatures T < T c we see that there is a sharp rise 
in voltage for the apparent critical current I C (T), which 
decreases with T, vanishing at T c . In Fig.l we plot the 
I C (T) line obtained with a voltage criterion of V < 10~ 4 . 
For currents below I C {T) there is a pinned vortex lattice 
(PVL) which is the same as the T = ground state, 
with delta-like Bragg peaks. On the other hand, for cur- 
rents I > I C (T) there is a moving vortex lattice (MVL), 



which has anisotropic Bragg peaks in the structure factor 
5(k) as shown in Fig. 2(a). There are two features in the 
anisotropy of 5(k) : (i) The height of the peaks decreases 
in the direction of vortex motion (i.e. perpendicular to 
the current drive) . (ii) The width of the peaks increases 
in the direction perpendicular to vortex motion. This 
means that thermal broadening is less effective in the di- 
rection of motion. We have also studied the behavior of 
the Bragg peaks of the PVL and the MVL for three dif- 
ferent lattice sizes (L = 50, 100, 150) and different values 
of I and T. In the PVL the peak height is independent of 
the system size, as expected for a pinned lattice. On the 
other hand, for the MVL, the peak height decreases with 
system size, with a power law behavior 5(G) ~ L~ r, G{ I >'T) 
and < rjc(I,T) < 2. In Fig. 2(b) we show a finite size 
analysis for two Bragg peaks S(Gi) and 5(d) at a par- 
ticular point in (I,T) with a power law fit. This is the 
expected behavior for a floating solid in two dimensions 
. In general we see that r\Q increases with T for a given 
current. When r] G (I,T) > 2 the MVL melts into a liq- 
uid. The anisotropic structure of the MVL of Fig. 2(a) 
is similar to the behavior predicted for a moving Bragg 
glass [0. However in our case there is no random pin- 
ning, but periodic pinning. The only source of dynamic 
randomness are thermal fluctuations. 

We now study in more detail the different transitions 
by fixing a given value of the current / and slowly chang- 
ing temperature T with small increases of AT = 0.0005. 
In this way we obtain the phase diagram shown in Fig. 1. 
We have also cross-checked these results with the IV 
curves at fixed T. There are two cases of interest: (i) 
low currents I < i c (0) and (ii) large currents I > I c (0). 

(i)Low currents. In Fig. 3 we show the behavior for 
/ = 0.03 < I c (0). At low temperatures the voltage is 
almost zero since the VL is pinned. When increasing T 
there is a sharp rise of the voltage (of two orders of mag- 
nitude) at a depinning temperature T P (I), as shown in 
Fig. 3 (a). At this temperature the VL starts to move since 
the driving current is higher than the critical current. 
Therefore, this corresponds to a transition from a pinned 
VL to a moving VL. We find that the I C (T) line obtained 
from the IV curves at fixed T coincides with the T P (I) 
line obtained from the V — T curves. We have also calcu- 
lated at the same time the structure factor S(k, T) and 
the transverse helicity modulus T^T). In Fig. 3(b) we 
show the behavior of two Bragg peaks S(G\) and S^G^). 
For T < T P (I) we see that 5(d) = 5(d) since there 
is a pinned VL with isotropic structure factor. Above 
T P (I) we find that 5(d) ^ 5(d)- This shows the fact 
that there is a MVL with anisotropic Bragg peaks. These 
peaks vanish at a higher temperature Tm(I) in a continu- 
ous and smooth transition, corresponding to a melting of 
the MVL. Above Tm(I) all Bragg peaks vanish and there 
is a vortex liquid for T > Tm- In the inset of Fig. 3(b) 
we show 5(d) for two different sizes L = 50,100 [we 
find similar finite size effects for 5(d)]- We see that for 
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T < T p (I) the S(Gi) is size independent since the VL is 
pinned, while for T p {I) < T < Tm{I) there is a power- 
law size dependence as expected for a floating solid, see 
Fig. 2(b). The temperature Tm{I) where S{G\) vanishes 
is size independent. On the other hand, the T X (T) has 
a significant drop at T P (I), however it remains finite but 
with large fluctuations in the MVL phase, T p < T < Tm- 
This suggests that in the MVL there is superconducting 
coherence in the direction transverse to the current. 

(ii) Large currents. When the VL is driven with a large 
current / > I c (0) there is a moving VL with anisotropic 
Bragg peaks even at T = 0. In Fig. 4 we show our 
results for / = 0.16 > / c (0). The structure factor is 
always anisotropic as can be seen in Fig. 4(b) where 
S{G\) ^ S{G2)- We find that the Bragg peaks vanish 
at a temperature Tm{I), which is size independent, as 
shown in the inset of Fig. 4(b) for S(G\) [similar be- 
havior is found for S(G2)} ■ We have investigated the 
possibility of a transverse critical current I c , tr f|Jl(J- At 
T = a finite 7 Cjtr is expected due to commensurabil- 
ity effects 10 7j. We have applied a transverse current 
hr = Ix in the x-direction (in addition to the applied 



bias, / 



0.16) and we have calculated the trans- 



verse voltage response Vt r . We find that for finite low 
temperatures Vt r is negligible small within our numerical 
accuracy, until there is a sharp increase at an "effective" 
transverse critical current 7 Cjtr (T). We find that I c ,tr(T) 
tends to vanish at a temperature T tr . An interesting way 
of studying this phenomenon is to apply a small current 
Itr and vary T. In Fig. 4(a) we study the onset of the 
transverse depinning transition: we apply a small current 
Itr = 0.01, and we show the transverse resistive response 
Vtr /Itr vs T. We see that Vt r is vanishingly small at low T 
and it rises at T tr . This transition temperature is clearly 
below Tm as we can see in Fig. 4. We have obtained 
T tr (I) for two driving currents / > / c (0) as we show in 
Fig. 1. It seems reasonable that this transverse depinning 
line will exist all along this region of the phase diagram. 
For the intermediate temperatures T tr < T < Tm we find 
that there is always an ordered vortex array but the ori- 
entation and structure of the MVL depends on the initial 
conditions. Moreover, finite size analysis shows that for 
T < T tr the MVL structure factor is weakly size depen- 
dent with t]g(T) ~ 0.05, while after T tr , the exponent rjo 
has a steep increase to values t]g(T) ^0.5. The helicity 
modulus is shown in Fig. 4(c). We see that for T < T tr 
the MVL has transverse superconducting coherence with 
a well-defined T^. On the other hand, T^, shows strong 
fluctuations in the region T tr < T < Tm- 

In conclusion, we have studied the current-temperature 
phase diagram of a vortex lattice driven in a periodic 
Josephson junction array. We find that for low currents 
the "longitudinal" depinning transition of the VL and 
the later melting of the moving VL are different with 
Tm(I) > Tp(I)- For large currents we find a very anal- 
ogous behavior for the "transverse" response of a fastly 



moving VL (compare Fig. 3 with Fig. 4). In this case it 
is possible to define a transverse depinning transition at 
a temperature T tr , and a later melting transition of the 
moving VL at Tm(I). This transverse depinning transi- 
tion could be easily studied in controlled experiments in 
Josephson junction arrays with transport measurements. 
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FIG. 1. I -T Phase diagram for / = 1/25. T M {I) line 
is obtained from T x vs. T curves (□) and from 5(G) vs. 
T curves (•). T P (I) line is obtained from IV curves (*), from 
5(G) vs. T curves (•) and from (V y ) vs. T curves (A). T tr (I) 
curve is obtained from (Vt r ) vs. T curves (o). Temperature 
is normalized by Io^o/^kB- 



FIG. 2. (a) Intensity plot of the structure factor 5(G) for 
a moving vortex lattice at I = 0.06, T = 0.02. (b) Finite 
size analysis and power law fit of 5(G) ~ L~ VG ; we obtain 
j? Gl = 0.25 ± 0.06 (*) and t)g 2 = 0.14 ± 0.05 (A). 



FIG. 3. For I < J c (0), I = 0.03: (a) Dissipation (V v )/I 
vs. T. (b) Structure factor at two lattice vectors, 5(Gi) (★) 
and 5(G 2 ) (o) vs. T. Inset: Size effect in 5(G). (c) Helicity 
modulus T x vs. T. 



FIG. 4. For / > J c (0) , / = 0.16: (a) Transverse dissipation 
(Vtr)/Itr vs. T. (b) 5(Gi) (*) and 5(G 2 ) (o) vs. T. Inset: 
Size effect in 5(G). (c) T* vs. T. 
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